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Abstract: We present ansatze that reduce the mass-deformed ABJM model to 
gauged Abelian scalar theories, using the fuzzy sphere matrices G a . One such reduc- 
tion gives a Toda system, for which we find a new type of nonabelian vortex. Another 
gives the standard Abelian-Higgs model, thereby allowing us to embed all the usual 
(multi-)vortex solutions of the latter into the ABJM model. By turning off the mass 
deformation at the level of the reduced model, we can also continuously deform to 
the massive 4 theory in the massless ABJM case. In this way we can embed the 
Landau- Ginzburg model into the AdS/CFT correspondence as a consistent trunca- 
tion of ABJM. In this context, the mass deformation parameter /x and a field VEV 
((f)) act as g and g c respectively, leading to a well-motivated AdS/CMT construction 
from string theory. To further this particular point, we propose a simple model for 
the condensed matter field theory that leads to an approximate description for the 
ABJM abelianization. Finally, we also find some BPS solutions to the mass-deformed 
ABJM model with a spacetime interpretation as an M2-brane ending on a spherical 
M5-brane. 
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1. Introduction 



Since its beginnings in 1997, the AdS/CFT correspondence [1] has found application 
in a variety of phenomena; not only in quantum gravity but also, increasingly in 
fields as diverse as low energy QCD and condensed matter. Its original formulation 
described four dimensional M = 4 Supersymmetric Yang-Mills theory with an SU (N) 
gauge group in the large N limit from the perspective of a dual gravitational theory 
on AdSa x S 5 . As a toy model for the exploration of four dimensional QCD at strong 
coupling, M = 4 SYM, with its large supersymmetry, conformal invariance, and 
large number of colors N (that ensures that the dual is just a gravitational theory 
and not a full string theory) is nearly ideal. By modifying this simple set-up, in 
particular by breaking supersymmetry and conformal invariance, a lot was learned 
about QCD itself as, for instance, in the Sakai-Sugimoto model [2]. A crucial part 
of this development is that the physics of gauge theories at finite temperature shows 
remarkable universality, which has translated into applications of M = 4 SYM to the 
high temperature plasmas at RHIC and the ALICE experiment at the LHC (see for 
example [3] for an extensive review and references). 

With the discovery of the pp-wave/BMN correspondence in 2002 [4] came the real- 
ization of the importance of operators with large R-charge to a full string theory (and 
not just supergravity) description of the dual. This, in turn, led to the description 
of spin chains from string theory [5] and, more generally, to an understanding of 
the integrable structures on both sides of the correspondence. In a sense, this was 
the precursor to the application of the gauge/gravity duality to condensed matter 
physics. More recently, another important example of AdS / CFT appeared, the so- 
called ABJM model; a three dimensional M = 6 supersymmetric Chern-Simons gauge 
theory with the gauge group SU(N) x SU(N), dual to string theory on AdS A x CP 3 . 
This model can be considered as a prototype for strongly coupled theories in three 
dimesions, in particular for planar condensed matter systems. For instance, in [6, 7] 
it was used to study the relativistically invariant quantum critical phase and com- 
pressible Fermi surfaces, respectively. These applications of the AdS/CFT corre- 
spondence to condensed matter hinge on the idea that, if physics in AdS is always 
holographic, then we can consider simple theories in AdS, which should be dual to 
some strongly coupled conformal field theories (see e.g. [8, 9] for a review). In an 
overwhelming majority of cases considered, the argument for applying the AdS/CFT 
duality (and trusting the answers it provided) was universality. In other words, the 
variety of theories usually considered contain a small subset of abelian operators dual 
to a small number of fields in AdS, usually a gauge field, some scalars and perhaps 
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some fermions. On the other hand, the relevant condensed matter models one usually 
wants to describe is usually abelian to begin with. It is not entirely clear then why 
we can either: i) focus on a small subset of abelian operators of a large N system; or 
ii) consider an abelian analog of the large N system, which would not have a gravity 
dual. 

A better motivated scenario for such an "AdS / CM T " correspondence would be if, in 
a large N field theory with a gravity dual, we could identify a consistent truncation 
of the (in general, nonabelian) field theory to an abelian subset corresponding to 
the collective dynamics of a large number of fields, and the resulting abelian theory 
would be a relevant condensed matter model. It is toward this end that we explore 
possible abelian reductions of the ABJM model in this article. Our strategy will 
be to look to the matrices G a that characterize the "fuzzy funnel" BPS state of 
pure ABJM and the "fuzzy sphere" ground state of the massive deformation of 
ABJM (mABJM) since they correspond to a collective motion of 0(N) out of 0(N 2 ) 
degrees of freedom. They will play a central role in our abelianization ansatz. We 
will then argue that this ansatz furnishes a consistent truncation of mABJM and 
can be used to identify further (phenomenologically) interesting abelianizations. We 
then show how these find application in condensed matter physics and, finally, we 
will explore some BPS solutions suggested by the abelian ansatze together with their 
spacetime interpretation. The main ideas about the abelianization and application 
to AdS/CMT were outlined in the letter [10], and here we present the full details. 

The rest of this paper is organized as follows. In section 2 we explore general abelian- 
ization ansatze involving G a , and identify two important cases of further consistent 
truncations for this model. In section 3 we study one of them, which, for BPS 
states, leads to a Toda system that possesses vortex-type solutions with topological 
charge and finite energy, but with \<f>\ — > at both r — > and r — > oo, which we 
describe numerically. In section 4 we describe a second case, more relevant for the 
AdS/CMT motivation above and find a reduction that, depending on certain param- 
eters, gives us either an abelian-Higgs model, or a 4 (relativistic Landau-Ginzburg) 
theory. In section 5 we study the relevance of this reduction for condensed matter 
and AdS/CMT and sketch a simple condensed matter model that reproduces the 
general features of abelianization. In section 6 we study some BPS solutions sug- 
gested by the abelianizations. Finally, in section 7, we provide a possible spacetime 
interpretation for these solutions in terms of M2-branes on a background spacetime. 

2. ABJM, massive ABJM and their Truncations 

The ABJM model [11] is obtained as the IR limit of the theory of N coincident 
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M2-branes moving in IR 2,1 x C 4 /Z fc . It is a M — 6 supersymmetric U(N) x U(N) 
Chern-Simons gauge theory at level (k,—k), with bifundamental scalars C 1 and 
fermions ipi, 1 = 1, ...,4 in the fundamental of the SU(4)r symmetry group. The 
gauge fields are denoted by and A^. Its action is given by 

S = Jd 3 x (j^e^ Tr (a^ v A x + jA,A u A x - A,d u A x - |i M i„i A 

- Tr D^CjD^C 1 - i Tr^'V^z + ^ Tr (c I C\C J C\C K C\ c 

+ c\c l c\c J c ] K c K + ac i c ] j c k c\c j c ] k - §c I c\c J c\c K tf K c K 

+ ^ Tr (cJcVVj - ^tfCtyj - 2C}C VVj + 2^ Jt C 7 C^j 



+ e™<^jCkVfc - e/m^^V 1 ) j (2.1) 

where the gauge-covariant derivative is 

D.C 1 = d.C 1 + xA^C 1 - iCjA^. (2.2) 

The action has a SU(4) x C/(l) R-symmetry associated with the M = 6 supersym- 
metries. 

There is a maximally supersymmetric (i.e., preserving all M = 6) massive deforma- 
tion of the model with a parameter /i [12, 13], which breaks the R-symmetry down 
to 377(2) x SU(2) x U(l) A X C/(1) B X Z 2 by splitting the scalars as 

C" = (Q Q ,iT); a = 1,2 (2.3) 

The Z2 action swaps the matter fields Q a and R a , while the SU(2) factors act 
individually on the doublets Q a and R a respectively and U(1)a symmetry rotates 
Q a with a phase +1 and R a with a phase —1. The mass deformation, besides giving 
a mass to the fermions, changes the potential of the theory. The bosonic part of the 
deformed action can be written as 

k ( 2i 2i * * * 

^Bosonic = |^ A Tr \A„d v A x + —A fl A u A x - A^d v A x - —A /1 A U A X 

- Tr \D^Q a \ 2 - Tr \D^R a \ 2 - V (2.4) 

where the potential is 

V = Ti (\M a \ 2 + \N a \ 2 ) , (2.5) 

and where 



M a = fiQ a + y (2Ql a g^ + R^R^Q a - Q a R^ + 2Q^R^R a - 2R a R^Q^ , 
AT" = + y (^R^R^R® + Q^iT - iTQ^ + 2R^Q\Q a - 2Q a Q) 3 R^ . 
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(2.6) 

The equations of motion of the bosonic Lagrangian (2.4) are 

dV dV 



dQV " dRi 

2vr , A 2tt 



where the field strength = d il A v — d v A^ + ([A^, A u ] and the two gauge currents 
J M and J M , given by 

J» = i (Q a {D»Q a y - (D»Q a )Q\ + R a (D^R a y - (D"R a )Rl) , 

(2.8) 

> = -i (Ql(D»Q a - (D^Q a ) ] Q a + Rl(D»R a - (D^R^R 01 ) , 

are covariantly conserved so that V M J M = V M J^ 1 = 0. In addition, there are two 
abelian currents j M and j M corresponding to the global U(1)a and U(1)b invariances, 
given by 

f = i Tr (Q a (D»Q a y - (D^Q a )Ql + R a (D t *R a ) ] - (D»R a )Rl) , 

(2.9) 

f = - tTr {Ql(D»Q a ) - (D»Q a yQ a + R ] a (D»R a ) - (D^R a )^R a ) , 

which are ordinarily conserved i.e. d^j^ = = 0. By choosing the gauge A = 
A = 0, the energy density (Hamiltonian) is given by 1 

H = Tr [(D Q a Y(D Q a ) + (AQ a ) t (AQ a ) 

+(D R a y (D R a ) + (DiR a )^ (DiR a ) + V] . (2.10) 

Since this is a Chern-Simons theory, the equations of motion must be supplemented 
with the Gauss law constraints 

F 12 = —J° = — (Q a (D°Q a Y - (D°Q a )Ql + R a (D°R a ^ - (D°R a )Rl) , 

(2.11) 

A /7T7 a Z.1T1 

F 12 = — J° = -— (Ql(^°Q Q ) - (D°Q a yQ Q + Rl(D°R a ) - {D Q R^R a ) . 
k k 

The gauge choice is not as restrictive as it would seem. Choosing (as we do below 
for our abelianization) A and A different from zero produces an extra term in 
the Hamiltonian of the form € fJ,uX Tr[A fl A u A\ — A^A^Ax]. In the abelian case this 
vanishes anyway since it is proportional to and there are only two 



1 Note that the terms A1A2 — A2A1 cancel from pq — L, and the rest of the CS term involve Aq 
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aft's. So in the abelian case, the Hamiltonian is the same even away from the gauge 
A = A = 0. The mass deformed theory has ground states of the fuzzy sphere type 
given by 

R a = cG a ; Q a = and Q\ = cG a ; R a = (2.12) 
where c = and the matrices G a , a = 1, 2, satisfy the equations [12, 13] 

G a = G a G\G p - G p G\G a . (2.13) 
It was shown in [14, 15] that this solution corresponds to a fuzzy 2-sphere. 

An explicit solution of these equations is given by 

(G l ) m)n = Vm - 1 5 m , n , 
(G 2 ) min = y /(N - m) S m+1>n , 

{G\) m ,n = Vm - 1 S m)n , 

{G\) m>n = V(N - n) S n+hm . (2.14) 

Clearly, these matrices satisfy G 1 = G\ also. In the case of the pure ABJM, there is 
a BPS solution of the fuzzy funnel type with c replaced by 



instead, where s is one of the two spatial coordinates of the ABJM model. The 
matrices G a are bifundamental under U(N) x U(N), therefore G 1 G\ and G 2 G\ are 
in the adjoint of the first U (iV), and G^G 1 and G\G 2 are in the adjoint of the second. 



2.1 An Abelianization Ansatz 



Given all these properties of the G a matrices, it is reasonable to choose the following 
abelianization ansatz 

A, = afG'G\ + a^G 2 G\ , 
A, = afG\& + flWGjG 2 , 
Q a = KG a , 

R a = Xa G\ (2.16) 

with no summation over a in the ansatz for Q a , R a ; a]L and real- valued vector 
fields and <p a , \a complex- valued scalar fields. 

Since G 1 G\ commutes with G 2 G\ and G\G l commutes with G\G 2 , the gauge fields 
aft are abelian and the field strengths decompose as 

F, v = d,A u - d w A fM + i[A„ A v ] = f^G'Gl + f^G 2 G\ , 

(2.17) 

F, v = d,A u - d v A, + z[4, A„] = fjSGlG 1 + f$G\G 2 , 
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with the abelian field strengths fj$ = d^ai — d v a$ . 
With this ansatz, the Chern-Simons term becomes 

while the covariant derivatives D^Q a and D^R give rise to 



(2.19) 



(2.20) 



and the values for M a , N a are given by 

M 1 = f [0 x (c 2 + |0 2 | 2 - | X2 | 2 ) - IX.XM G l , 
M 2 = ^ [0 2 (c 2 + |0 X | 2 - | X i| 2 ) - 2xIX20i] G 2 , 

iV 1 = T [xi( - c 2 + |x 2 | 2 - |0 2 | 2 ) - 2fcfc X2 ] G 1 , 

iV 2 = y [ X2 ( " c 2 + |xi| 2 - I0i| 2 ) - 2^02X1] G 2 
where as before, c 2 = fik/ (2n). Substituting into the potential gives 

V = 2 -^N(N - 1) [(l^l 2 + |xi| 2 )(|x 2 | 2 - |0 2 | 2 - c 2 ) 2 

+ (|0 2 | 2 + |X2| 2 )(|xi| 2 -|0i| 2 -c 2 ) 2 (2.21) 

+ 4|0 1 | 2 |0 2 | 2 (|xi| 2 + |X 2 | 2 ) + 4|xi| 2 |X2| 2 (|0i| 2 + |0 2 | 2 )" 

Note that the interchange of x with <fi (which changes Q a with FT) is equivalent to 
a change in the sign of c 2 , i.e. either a change in the sign of //, or of k. Putting 
everything together then gives the final abelian effective action 



' d 3 x 



/• 



47T 

(2.22) 

with a rescaled potential U = 2V/N(N — 1). Since the effective theory derives from 
a Chern-Simons theory, the equations of motion need to be supplemented with the 
Gauss law constraints which, in our ansatz, reduce to 



fi? = y[^iW " (W! + XiP°Xi) f - (D° Xl )x\} , 



/i? = ^[0 2 ( J DV 2 ) t - P°0 2 )0 2 + X2( J D°X 2 ) t - (D° X2 )xl} , 



(2.23) 
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We see, however, that these are nothing but the equations of motion for the 

action (2.22). As we will need to work away from the = = gauge, we don't 
need to impose them. 

2.2 Consistent Truncations 

A key point to note about this abelianization ansatz is that it a consistent truncation 
of the original ABJM theory in the sense that, using the facts that M a oc G a , 
N a oc G«, D^{<j> a G a ) = (D„D^ a )G a and D^( Xa G a ) = (D^ X a)G a , the 
equations of motion that follow from the action (2.22), 



k_ 

^_ ,w\A2) _ 

4vr ^ " 



b 2 D x <p 2 - <p 2 D^<p 2 + x 2 D x X2 ~ X2D X X 2 



(2.24) 



and 



4tt 2 r 

ir[(lx 2 | 2 
(Ixil 2 



|2 2\2 



Y + 2(|0 2 | 2 + | X2 | 2 ) (|^| 2 + | Xl | 2 + c 2 ) + 4\fo\ 2 \ X 2\' 



4vr 2 



1 2 2\ 2 



fc 2 L 
4vr 2 r 

^[(lx 2 | 2 
r (lxil 2 



c 2 )^ + 2(|0 1 | 2 + | Xl | 2 ) (|0 2 | 2 + | X2 | 2 + c 2 ) + 4|0 1 | 2 | Xl ' 2 



/ |2 2"* 2 



4vr 2 



k 2 I 



) +2(|0 2 | 2 + | X2 | 2 )(|0 1 | 2 + | Xl | 2 -c 2 )+4|0 2 | 2 | X2 | 2 
x | 2 - c 2 ) 2 + 2(|0 1 | 2 + | Xl | 2 ) (|0 2 | 2 + | X2 | 2 - c 2 ) + 4|0 1 | 2 | Xl | 2 



Xi 



X2 , 

(2.25) 



satisfy the higher original ABJM equations of motion (2.7) and Gauss constraints 
(2.12). 

Since Tt[G 1 G\] = Tr[G\Gi] = Tr[G 2 G 2 ] = Tr[G^G 2 ] = N(N - l)/2, the energy 
density (Hamiltonian) is 

H = N ^ N ~ V + \D oXt .\ 2 + |I>„&| 2 + \D aXi \ 2 } + V (2.26) 



where a, b = 1,2. Note also that away from the gauge A = A = (which imply 
that Oq = 0), we would, in principle, have a term cubic in the gauge fields in the 
Hamiltonian. This however vanishes in the abelian case, so the above result is correct 
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in general. This abelianization, with its four complex scalar fields, is rather general. 
We will study further reductions of it involving only two scalars. Looking at the 
scalar equations of motion above we see that putting any two of the scalars to zero 
is again a consistent truncation. 



A trivial choice turns out to be \2 = 02 = (or equivalently Xi — 4>i = 0)> 
since in that case, the potential reduces to a simple mass term, 

V=^m 2 + \xi\ 2 ) (2.27) 

(2) 

while at the same time, the only ct/t dependence remains in the Chern-Simons 
term, ~ J ea^f^\ so its equation of motion is ff$ = 0, which means 
is also trivial (pure gauge). So we remain with two massive complex scalar 
fields coupled to one trivial gauge field (pure gauge, with no kinetic term), an 
uninteresting model. 

A much more interesting choice is 0i = 02 = 0, which will turn out to lead 
(with some modifications) to the Abelian-Higgs model. Since we will study 
this separately and extensively in section 4, we will not discuss it further here. 

Finally, setting \i — 02 = 0, and renaming \2 to 2 for simplicity, we get 



s = - N < N -V U X 



^ x {tf ] fS! + ^fS) + \DM 2 + u(\U N) 



(2.28) 

V = 2 ~^N{N - 1)[N 2 (|0 2 | 2 - c 2 ) 2 + |fc| a (|fc| 2 + c 2 ) 2 ] , 
and energy density (Hamiltonian) 

H = N(N ~ l \ \D^ 2 + lA^I 2 ] +V. (2.29) 



We should note that, until now, we have worked only with the massive ABJM model, 
but that we can analyze the massless (or pure) ABJM model in a straightforward 
way by setting c = 0. Since c appears only in the potential, we can check that the 
model with potential (2.21) is symmetric under interchange of <pi -H- x%- F° r the 
model with xi — 02 = above, for example, we obtain a purely quartic potential, 

2tt 2 

V = ——N(N - 1)|0 1 | 2 |0 2 | 2 (|0 1 | 2 + |0 2 | 2 ) . (2.30) 
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3. New vortex solutions for a Toda system 



We now study BPS solutions of the effective model (2.28). Before doing so, it is 
worth taking a step back, and considering the more general case of the Q 2 = R 1 = 
reduction, with only Q l = Q and R 2 = R nonzero, but without the abelianization 
ansatz. There we can 'complete squares' in the Hamiltonian density and write it, in 
complete analogy to the usual Abelian-Higgs model, as 



H = Tr \D Q - iMf + Tr \D R + iN\ 2 + Tr \D^Q\ 2 + Tr \D + R\ 2 

+ ie ab d a Tr (Q\D b Q) - R\D h R)) + /ij , (3.1) 



where fijo = fik/ (2-n) Tr (-Fi 2 ) and D± = D\ ± iD 2 . Just as in the Abelian-Higgs 
model, the term on the second line (with e ab ) is zero on the configurations of interest, 
since 



/ d 2 xe ab d a (^D l 
Jv 



] D a (f)) dx a 



(3.2) 



and D a (p — > at r — > oo for <p = Q, R in order to have finite energy configurations. 
Moreover, the perfect squares on the first line are minimized by the BPS equations 



D_Q = 0; D + R = 0; D Q = iM\ D R 



-iN . 



(3.3) 



which leaves just the topological term, fijo- The BPS equations together with the 
Gauss law constraints are 



D_Q = 0, 
D+R = , 

D Q = ifiQ 



DqR = ifiR 



F\2 
F\2 



2m 



2m 



47T/X 



QR ] R - RR ] Q 



RQ^Q - QQ^R 



7T 



(3.4) 



AlTfl 



8tt s 



QR ] RQ ] - RQ ] QR ] 
R^QQ^R - Q^RR^Q 



where, in the Gauss law constraints, we have already substituted the BPS equations 
for D Q, D R. These equations are more general, and can be used in principle to find 
nonabelian BPS solutions. In practice, they are still too difficult to solve analytically 
so from now on we will go back to the abelian case Q = (ftiG 1 , R = faG 2 . There, the 
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BPS equations for D Q and D R become 

2vri 



(d - ia^)^ 



k 



fj,k 
2^ 



(3.5) 



{do 



i a 



(2)> 



2m 



fik 
2^ 



For static configurations (for which do<pi = 0) these equations can be solved for ag 



,(2) 



2tt 



2^ 



(3.6) 



2tt 

T 



In other words, the ag are completely specified by the scalar fields <pi and spatial 
components of the abelian gauge fields da , a = 1,2. Consequently, the (temporal) 
gauge (Zq = would be inconsistent with the BPS equations. This is different from 
the Abelian-Higgs model, where the one can set both do = and 8q — 0, reducing 
the system to a two dimensional one (for the spatial components). Here, this would 
be inconsistent with the Gauss law constraint which, for a Chern-Simons gauge field, 
relates F i2 to terms with D Q and D R so that, if F i2 is nonzero, so too is D Q Q and 
D R. Finally, the Gauss law constraints in the BPS case reduce to 



f (i) 

J 12 



f (2) 
J 12 



57T 



k 2 



8tt^ 
k 2 



>i| 2 + ^ 
In 



fik 
2^ 



(3.7) 



In order to facilitate the rest of the analysis of the BPS system, it will prove useful 
to complexify the (xi, x 2 )— plane and write z = x l + ix 2 . As usual, this induces a 
complexification of the derivatives as well as the gauge fields as, 



8=^(4 



id 2 ) , a (i) = 



.(*) 



lay) 



together with their complex conjugates. This, in turn, implies that f 
Ba^}, so that the BPS equations D_Q = D + R = become simply 

90i - ia {1) 0i = 0, d(p 2 - ia {2) <p2 = . 



(<) 

12 



-2i[9oW- 



(3. 



Equations (3.8) together with the Gauss law constraints constitute a complete set 
which, as we argue below, possess at least one simple set of finite energy, spatially 
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localized solutions of the vortex type i.e. isolated zeros of the (complex) scalar fields 
with nonvanishing winding number. The analysis follows the same general logic as 
for the Nielsen-Olesen vortex and we start by writing 

1 = |0i|e iei ; = (3.9) 
and then (3.8) become (after taking derivatives and making the combinations ) 

2dd\n\<p 1 \ 2 - 2i(8d - dd)6 l = ^A In |0x| 2 + e ah d a d b 6 1 , 

(3.10) 

2dd In |0 2 | 2 + 2i(8d - dd)6 2 = ^Aln |0 2 | 2 - e ab d a d b 9 2 . 

But since, if a is the polar angle in the 1,2 plane, e ab d a dbCt = 2ir5 2 (x) (as can be 
checked by integrating over a circle of vanishingly small radius), we may take the 
ansatz 

6 1 = -N 1 a, 9 2 = N 2 a, (3.11) 

which leads to 

^A In |0! | 2 -2ttN 1 5 2 (x) , 

(3.12) 

^Aln|0 2 | 2 -2vriV 2 «5 2 (x). 

Finally, on substituting the Gauss law constraints, we obtain a continuous Toda 
system with delta functions sources, 

A In |0i 
A In |0 2 

whose solutions we proceed to analyze. 

3.1 Asymptotic Analysis of the Toda System 

As in the case of the (much simpler) Nielsen-Olesen vortex, the Toda system of 
equations does not, as far as we are aware, exhibit any closed form analytic solution. 
Consequently, here too must we resort to topological, asymptotic and numerical 
analyses to tease out finite energy solutions from it. The argument, fortunately, goes 



f(i) _ 

f(2) _ 
J 12 ~ 



f(l) _ 

J12 ~ 

f(2) _ 

J 12 ~ 



~k 
Air 

T 



fik\ 



2n J 



+ 4nN 1 5 2 (x) 



(3.13) 



Ztt 
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through in much the same way as for the Neilsen-Olesen case: the topological term 
in the energy 

f , , f , , ^ m /r, n ukN(N-l) [ , , , ri 9 \ 
/ dxdyfxjo = / dxdy — Tr(F 12 ) = / eta A/ (f 12 + / 12 ) , 

is quantized as usual, since for an abelian gauge field 



!/"_.. 1 /" . „ 1 



27T ./nt2 27T ./c 27T ,/ 



F 12 dxdy = — (b A a dl — — I A a da , 



2tt 



and D a (j) — > at r — >• oo, with = l^le* 6 *, = — za M and c^ln |0| — >■ 0, means 
that d a 0°° — = and consequently 

F 12 dxdy = 2-kN . 

The corresponding statement for our system (3.8) is that, at infinity 

« floo (i)oo 

which gives the energy of the BPS state as 

N(N - 1) 

E(N 1} N 2 ) = yk-^ ^(JVl + N 2 ) . (3.14) 

We are now in a position to look at the Toda equations in the asymptotic regions. 
To obtain the r — > behaviour, we integrate each of them over a very small disk of 
radius R — > 0, and find that 

dxdy V • Vln = 2^^ , z = 1, 2 , 



which, after using Stokes' theorem, gives 

i?— ln|0j|| r=i? = iVj. 

This expression is easily integrated to show that as r — y each of the scalars exhibits 
the power law behaviour, 

\<f>i\ ~ Air"' . (3.15) 

This is in accordance with the usual argument says that the only possibility for the 
vortices with <fi = |0|e i7VQ is to have \<f>\ — y at r — > in order that the phase is well 
defined at r = 0. In fact we can do better and refine the conditions at r — y by 
using the equations of motion away from r = 0. Taking as an ansatz for the scalars 

|&| a = A l r 2N ^l + B 1 r^), 

(3.16) 

|0 2 | 2 = A 2 r 2N *(l + B 2 r«), 
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and substituting into the Toda equations, we find 



q = 2N 1 + 2, p = 2N : 
kp 



* = ^. (3.17) 



B 2 = _ 8 !^ j4l 

kq 2 

The constants are only determined from the full numerical solution. 

At r — > oo, we can first check that neither a constant, nor a decaying exponential, 
nor a power law that blows up, ~ r p works for either of the two fields. This leaves 
a decaying power law as the only plausible behaviour for either of the two scalar 
fields. In order to use the equations above, we must consider also the first subleading 
terms, i.e. we substitute 



M 2 = -(i + - 

|r 1 r n \ r q 



(3.18) 



into the equations above, to find 



m = q + 2, n = p + 2 , 

87T/X 

p 2 k 



B 1 = -(^)A 2 , (3.19) 



Since p,q G N*, m,n — 3,4,5,... Again, the constants Ai,A 2 , as well as m, n are 
determined from the full numerical solutions. 



3.2 Numerical Analysis of the Toda System 

To determine the various parameters of the vortex-like solutions described above, we 
need to solve the Toda system numerically. As in the asymptotic analysis above, our 
numerical solution follows the general logic of the Abelian-Higgs model. Specifically, 
we will use a modified two-parameter shooting method to numerically solve the two- 
point boundary value problem described by the coupled Toda equations. To facilitate 
the implementation of the shooting algorithm, we first rewrite the equations as a 
four-dimensional (non-autonomous) dynamical system. To this end, we first non- 
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dimensionalize the system by rescaling our variables and denning 



2 



f^-AfoW (3.20) 

with r taken to be the radial coordinate on the plane. Substituting into the system 
(3.14) and assuming that the solitons that we are looking for are rotationally sym- 
metric on the plane (so that the two-dimensional Laplacian is A = (r-r-)), we 
find that eqs.(3.14) reduce to 

1 d (Rdf 



RdR \ fdRj ~~ 9 ^ ^' 

(3.21) 

1 d ( R dg 



RdR\gdRj ~^ 9 + l > 

Finally, we reduce the order of the system by one by making the additional defini- 
tions h = and j = j% so that (denoting by a ' derivatives with respect to the 
dimensionless radial variable R) 

f = h, 

h'=%---gf{f-l), 
J r 

9 = J, (3-22) 

f = -- l -gf(g + i). 

9 r 

Before directly integrating this four-dimensional non-autonomous dynamical system, 
it will be instructive to extract some qualitative information from it. There are two 
(physical) fixed points at (f,h,g,j) = (0,0,0,0) and (1,0,0,0). A linearization of 
the system near the former, shows that the origin is a saddle. Solutions of the kind 
that carry nonvanishing winding number and conform to the asymptotic boundary 
conditions /(0) = g(0) = and /(oo) = g(oo) = correspond to homoclinic orbits 2 
that begin and end at (0,0,0,0) and that encircle the fixed point at (1,0,0,0) (see 
Fig.l). 

Our numerical integration of the system is based on a two-parameter shooting algo- 
rithm that converts the nonlinear dynamical system above into a nonlinear parame- 
ter estimation problem. The parameters in question are precisely the undetermined 

2 This should be compared to the standard ANO vortices of the Abelian-Higgs model which cor- 
respond to heteroclinic orbits interpolating between the two fixed points of the associated dynamical 
system. 
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Figure 1: The (/, /')— subspace of the full phase space of the Toda system 



constants A% and A 2 above and these are chosen at R = so that the constraint 
/(oo) = g(oo) = is met. In practice, the constraints at R = oo are a problem, but 
our asymptotic anaylsis above can be extended to show that solutions at R ~ 10 
are quite safely in the far field for both / and g. Some results of our numerical 
integration are presented in Figures 2 and 3. 



Figure 2: The N\ = N 2 = 1 soliton profiles. Optimization of the shooting parameters 
yeild A 1 = 30.00, A 2 = 30.05 



Figure 3: The N\ = 2,N 2 = 1 soliton profiles. Optimization of the shooting parameters 
yeild A 1 = 50.00, A 2 = 100.00 





N± = 2, N*2 = 1 
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We also obtain from the numerics that the power law at infinity is \4>i\ 2 oc 1/r 3 , |0 2 | 2 oc 
1/r 2 , i.e. m = 3,n = 2. 

As a final point, we note that in the massless ABJM case, at fi = 0, these vortices 
vanish since their energy is proportional to fi. This agrees well with known facts 
about the solitonic spectrum of pure ABJM [16]. 



4. The Abelian-Higgs model from ABJM 

We now look to embed the Abelian-Higgs model in ABJM, as a truncation of our 
general abelianization ansatz. To find the truncation we look at the multi-vortex 
solution we found previously in [17] for the N = 2 case, i.e. U(2) x U(2) ABJM. 
There, not only was the ansatz written in a manner similar to the multi-vortices of 
the conventional Abelian-Higgs model, but the action on the moduli space of vortices 
was also found to be the same. In retrospect, this was really a telling signal that 
we were actually embedding the Abelian-Higgs model into ABJM. For the reader 
unfamiliar with [17], we recall that the static multivortex solution there was given 
by 

C 2 = ^-* /2 ffo(s) (° J) = ^-*' 2 H (z)G 2 ; C s = C = 0, (4.1) 
4>=-W); ' P " ) : A, = A, <° ° 



where Hq(z) = Y\™ =1 {z — z{) is an arbitrary polynomial and the real function ip{z) is 
determined through the equation 

ddip = fi 2 (1 - e-^\H (z)\ 2 ) (4.2) 

with boundary conditions at |z| — > oo requiring ip — > log \Hq(z)\ 2 . As usual, Z{ with 
% — 1 . . . n, denotes the positions of the n vortices. Treating each of these position 
variables as (adiabatic) functions of time, zi(t), produces the first order solution 



c 1 


= o, 




— ^0 






4 1} 


= ( 2M 





J) -i{z 1 d l -£ i di)ijj 

#dgfy0\ i<i) = /0 o 
J ' z ' I ±z%diP ) ' 



(4.3) 
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on the moduli space of the vortices. On the other hand, when N = 2, our general 
abelianization ansatz gives 



= a% ) G 1 G\ + afG 2 G 



a? 



A„ = a^G\G l + afG\G 2 



a 







(2) ' 



(4.4) 



a 



(1) V aL 2) 



Comparing with the solution above (and also denoting now the first order solution 
for the abelian fields with a tilde to avoid confusion with the indices (1) and (2) on 
the a's) we find 



,( 2 ) 



4 2) = --(i^-!7^. 



dip, 4 1} = -<9<9</>, a« = a[, 2) = 



(4.5) 



Note that from (4.5), by taking complex conjugate and then sums and differences, 
we get 



s(i) 



1 

2/x 



+ I £?i - i ( - z l d- { ) d- 



(4.6) 



2 

2/j L 



)d 1 -i[z i d~ i + m )d: 



4>. 



This can be written more compactly, by using the fact that ip is a real-valued field, 

as 

(4.7) 



a„- - 



In view of the above solution, and assuming that the same relation to our abelian- 
ization holds at all N, we can now identify the truncation ansatz needed to obtain 
the abelian-Higgs model as 



<f ) 2 — 0, Xi = b = constant 



(4.8) 



which gives 



D 



fi<P2 = 

-ia^b 

A 1 



(4.9) 



D^X2 = (Dp - ia^')X2 



,(2)> 
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and the potential 



2tt 2 

y = _ ATfAT 1\[\U2n, |2 2,2 
k 

2tt 2 



2 iV(iV-l)[|6|^(| X2 |^-c^ + |x2r(H 2 



|2\2n 



k 2 



N(N - l)[\b\ 2 \ X2 \* + \ X 2\ 2 (-m 2 c 2 + \b\' + c 4 ) + c 4 |6| 2 ] 



(4.10) 



We can easily arrange for the coefficient of the \x\ 2 term to be negative, as is required 
for the mexican hat potential of the abelian-Higgs model, by choosing for instance 



|6| = |c| =>• jJL 



27r|6p 



(4.11) 



The action is then 
s _ N(N-l) 



a fx 



with the auxiliary field ofi . As usual it can be eliminated through its equation of 
motion 

k ,„.,..,., (4.13) 



2tt 



e „A a (l) f (2) 



(4.12) 



,(1) 



so that 



S 



N(N — 1) 



d x 



4tt|6| 2 

k 2 
8tt 2 |6| 2 



(4.14) 



which is nothing but the action of the abelian-Higgs model. 



Of course, we still need to check the consistency of the truncation, i.e. to check that 
the equations of motion of the full abelianization ansatz in section 2 are satisfied. 
We have fixed 0i, 02 to zero and xi f° b, so it is these three equations of motion that 
we need to check. As before, the choice <pi = <p2 = is a consistent truncation. The 

(2) 

equation for xi reduces, in the Lorentz gauge d^a^ = 0, and using (4.13), to 



dV 



(4.15) 



which is just the equation of motion we would obtain for the parameter b by varying 
in the abelian-Higgs action (4.14). We find this somewhat puzzling, since it means 
that the constant parameter \b\ has to be effectively treated like a field in the abelian- 
Higgs action, giving its own equation of motion. 

It remains now to check that our multivortex solution satisfies the condition (4.13), 
since it certainly matched our ansatz before we imposed the equation of motion for 
aft. The equations (4.13) reduce for the zeroth order solution and the first order 
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solution respectively, to 



-^e ij da {2) 



(4.16) 



^ = 



4n\b\* e dja ° 



(2) ' ~ 

We can check the first equation, since a\ = ^dip, which written in real components 
reads a!p = —^eijdjip. Then the zeroth order equation is satisfied if 



H = 



2ir\b\ 



k 



while the first order equation is satisfied when 



= -V%)a,-v> 



(4.17) 



(4.18) 



as expected. The appearance of the abelian-Higgs model above is somewhat non- 
standard but can easily be put into canonical form by appropriately normalizing the 



fields 



to obtain 



a 



(2) _ ^„(2) 

~ Nk a 
1 



X2 



X2 

N ' 



(4.19) 

where now = <9 M — igaffi and g = In terms of the canonical fields and 

coupling, the potential 



V 



fcl 4 + ^ + lx 2 | 2 iv 2 



Auk ,, l9 u 2 k 2 
+ \ b \ + ' 



2tt 



47T 2 |6I 2 



(4.20) 



As previously alluded to, the potential has a range of values of |6| for which it 
is spontaneously breaking (has negative mass squared). The central value of this 
domain is \b\ = c, and for this value of |6|, we find 



,„ ,o ukN 

M 2 - -^ 



21 



(4.21) 



Moreover, for this value of |6|, the equation of motion for \b\ 2 (the extra constraint 
on our abelian-Higgs model), becomes 



[\<h?-^ = \ 



A_ f (2)' " 
J u 



(4.22) 
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(2) 

On the other hand, equating the kinetic (Maxwell) term for f^J with the potential 
term, V, gives exactly the same equation. Further, taking the square root of this 
equation, and imposing that/ « = 0, we find 

^/r? = ±[|X2| - JW] , (4.23) 

which is part of the abelian-Higgs BPS condition. In other words, the extra condition 
is satisfied on BPS solutions of the abelian-Higgs model with f 0i = 0, and in particular 
for vortices. 

This completes our demonstration that (i) the abelian Higgs model can be obtained 
from the abelianization of the ABJM model and (ii) that both the classical (zeroth 
order) and the first order (in the moduli space approximation) multivortex solutions 
of the latter at N = 2, are encoded in this model. Since this is a bone fide embedding 
of the abelian-Higgs model, we say say more even. For instance, it is natural that we 
obtain the same fluctuation action for vortex scattering as in the abelian Higgs case. 
It also means that we can now immediately write down the multivortex solution 
at general N, with the guarantee that we will recover the same fluctuation action 
for vortex scattering as in the abelian Higgs case. To be concrete, the multivortex 
solution at general N in ABJM 

Ri = ^¥ G \ Ri = ^e-^ 2 H (z)G 2 , Q 1 = Q 2 = 0, 
A = -ddif>G 1 G\, A = -dd^G{G\ (4.24) 
A- z = l -B^G 2 Gl A- z = % -d^G\G\ 
produces an effective Lagrangian on the moduli space 

2 7T J [ 2 V f 1 J 

N(N-l) K _ i ^ k ^2_ km y^ KMz i_ zj ^.i_.^2 (J2rj) 



n 

-k^n + S ^—\z l \ 2 -k^q S ^K Q (2^,\z' - y )|-' - : 

i=l i>j 



with q ~ 1.71. To close this discussion on vortices of the abelian-Higgs model 
and their embedding into the ABJM model, we mention briefly that in the case of 
the massless ABJM model, with c = /x = 0, we obtain a non-symmetry breaking 
potential, 

2tt 2 

V = ~^N(N - 1) [|6| 2 | X2 | 4 + N 4 |X 2 | 2 ] , (4.26) 
which is just a massive gauged (f) 4 model. 
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5. Towards a string construction of AdS/CMT 



At this point, let's stop and consider what it is that we have achieved. Stripping 
away all the bells and whistles, essentially our truncation has produced a (2+1)- 
dimensional scalar field theory with potential 

V = 2 -^N(N - 1) [|fo| 2 |0| 4 + \<P\\(\b\ 2 " c 2 ) 2 " 2c>| 2 ) + c 4 |6| 2 ] . (5.1) 

It is not too difficult to see that it is just a Landau- Ginzburg model in which, at fixed 
\b\ 2 , c 2 oc /i acts as a coupling that takes us from a |0| 4 theory (the insulator phase) to 
an abelian-Higgs theory (the superconducting phase). In this sense, the parameters 
\b\ 2 and c 2 control the coupling g and critical coupling g c of the Landau-Ginzburg 
model. More precisely, we identify the combinations (|6| 2 — c 2 ) 2 as g and 2c 2 |6| 2 
as g c respectively. In this light, it makes sense then to think of this abelianization 
as a realization of the recently proposed AdS/CMT correspondence. To see why our 
construction is markedly different from any of its pre-cursors, we recall the general 
ideas involved. Usually, in an AdS/CMT construction, one assumes some theory in 
an AdS background, usually involving gravity, a gauge field A^, maybe a complex 
(charged) scalar and some fermions ipi. It is then argued that this theory should be 
dual to some large N conformal field theory with a global current J M dual to the gauge 
field A^, and some other operators (in principle) dual to the other fields. It is then 
argued that relevant physics in AdS corresponds to some behaviour of the operators 
in the field theory which simulates the relevant physics, like superconductivity [18] 
for example, to be studied. Sometimes the AdS theory is obtained as a consistent 
truncation of some known AdS / CFT duality (for which there is a heuristic derivation 
involving a decoupling limit of some brane constructions), so that the field theory 
contains a small subset of operators that could possibly give the desired physics [8, 9]. 

However, even in these cases, it is not obvious how to directly relate the set of 
operators in the given CFT to the condensed matter system of interest, and usually 
one has to invoke some sort of universality argument. In other words, if the physics 
of the selected set of operators in the large N CFT describes the correct physics for 
the condensed matter system, then perhaps the physics is general enough to appear 
in many different systems, and we can try to apply our seemingly unrelated field 
theory to the condensed matter system of interest. While we certainly appreciate 
the logic of this argument, we find it less than satisfactory for a number of reasons. 
Primary among these is that it is not at all clear why can we choose only a very small 
number of operators in the large N CFT and concentrate on their physics. Secondly, 
if we try to write down a gravity dual of an abelian theory having this small number 
of nontrivial operators, we would fail, since the absence of the large N would mean 
that we could not focus on the supergravity limit in the dual. 
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However, we can now do better. We have found a consistent truncation of the large 
N CFT, for which there is a well-defined duality, and not just a truncation of the 
gravity theory. That means that this set of fields is a well defined subset at the 
quantum leveP corresponding to the collective motion of the nonabelian fields in the 
large N case and involving 0(N) out of the 0(N 2 ) fields of ABJM, via the nontrivial 
matrices G a (which have O(N) nonzero elements). It is not just a simple restriction 
to AT = 1 of the ABJM model, which would imply losing the supergravity limit in 
the dual. Therefore this abelianization still maps to a purely gravitational theory, 
and not a full string theory as for generic abelian theories. 

In our case, there already exists a well defined gravity dual of the field theory. In the 
case of massless ABJM, that theory corresponds to M2-branes moving in the space 
M 2 ' 1 x C 4 /Zfc, and the gravity dual (i.e. the near-horizon limit of the backreacted 
background) is AdS± x CP 3 . In the case of the massive ABJM, the theory corresponds 
to M2-branes moving in a space defined in [17, 19] with the gravity dual described 
in [20, 17]. Of course, we still would need to understand to what the truncation 
t° (xi) = b and \2 ^ corresponds in this gravity dual in order to complete the 
picture, but we leave this for further work. 

Actually, as it turns out, the theory we obtain in the abelianization is also the relevant 
effective theory for a CMT construction. Indeed, as reviewed for instance in [21], 
starting from the Hubbard model for spinless bosons hopping on a lattice of sites % 
with short range repulsive interactions, 

H b = - W j2 ( b l b J + fyi) + 2 Yl n ^ - x ) ~ t* Ui ' ( 5 - 2 ) 

(ij) i i 

where = b\bi and w is the hopping matrix between nearest-neighbour sites, one 
obtains the relativistic Landau- Ginzburg theory 

S = J d'x^-m' + v'lV^ + ig-gcM' + uW 2 ) . (5.3) 

The effective field is obtained as follows. The ground state contains an equal num- 
ber of bosons at each site, with the creation operators a\ producing extra particles 
at each site, and creation operators h\ that produce extra "holes" at each site; "an- 
tiparticles" in the QFT picture. Then, as is usual in field theory, 4>i ~ cti^i + P%h\ is 
a discretized version of the complex field describing both particles and antiparticles, 
where aij, are wavefunctions for the modes. 

For g < g c we have an abelian-Higgs system, i.e. a superconducting phase, while 
for g > g c we have an insulator phase. The marginal case g = g c is a conformal 

3 We can consistently put the other fields to zero even at the quantum level. 
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field theory. The systems described by the above model also have a quantum critical 
phase which opens up at nonzero temperature for a T-dependent window around 
g = g c . This quantum critical phase is strongly coupled and very hard to describe 
using conventional condensed matter methods, which makes it an excellent choice 
for a holographic description. In [6] it was shown that by considering a gauge field in 
the gravity dual of ABJM and introducing a coupling for it to the Weyl curvature, 
7 / C a b c dF ab F cd one obtains a conductivity ct(oj) consistent with the quantum critical 
phase, and from which it was concluded that ABJM is a good primer for these 
systems, though the precise reason for the match was not obvious. 

While the bosonic Hubbard model leads, in the continuum limit to the action (5.3), 
the model itself is a drastic simplification, of a condensed matter system. The model 
has been used to describe the quantum critical phase of (bosonic) 87 Kb cold atoms 
on an optical lattice, but the description is believed to hold more generally for the 
quantum critical phase. For instance, high T c superconductors have a "strange metal" 
phase that is believed to be of the same quantum critical type. We can consider a solid 
with free electrons (fermions, perhaps several per atom) that could hop between fixed 
atoms, and unlike the simple Hubbard model, we also have in principle interactions 
that are not restricted to nearest neighbours. One could, for instance, generate 
bosons (pij (having the role of the bosons bi of the Hubbard model) by coupling 
fermions at two sites i and j. By an abuse of notation we will call by the same <ftij 
the field obtained by multiplying the corresponding "particle creation" operator with 
a wavefunction, and adding a corresponding "hole" part. 

In fact, we can sketch a simple model for the condensed matter system above that 
generates the same qualitative picture as the abelianization of the ABJM model. 
Consider spinless bosons 0^ generated by coupling fermions of opposite spins (Cooper 
pairs) at sites i and j with a maximum distance between sites \i — j\ < N, i.e. i/Jii/jj. 
The resulting faj can be described by a field with a,b = 1, ...N. Since we are 
in two spatial dimensions, every site has 0(N 2 ) neighbours a distance < N away. 
Now take the point i' at which the effective field, <pij, lives to be midpoint of the line 
between i and j, and a, b to correspond to sites j in the x and y directions away from 
i' (so that, if i' and j are fixed, so is i). Consider that the normalized wavefunctions 
for the field give probabilities for existence of the pairing as oc |<% 6 | 2 for a pair 
(ab). In this case, any transformation on <ftf must be a unitary transformation 
Uab,a'b> ingide u(N 2 ), up to an overall factor. In particular, any symmetry of the 
system must be of this type. The symmetry of the ABJM model is U(N) x U(N), 
and would correspond to U ab ' a ' b ' = fU aa 'V w '. 

Since the simplest type of condensed matter system is a rotationally invariant one, we 
should not have any angular dependence, and we should have = <pi{\/a 2 + b 2 ) = 
<p ba . It should be then possible to diagonalize this symmetric matrix, corresponding to 
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considering only the constant (rotationally invariant) m = modes for the "spherical 
harmonics" expansion e 2mm9 at fixed radius r = \/a 2 + b 2 . In this way, only O(N) 
modes, specifically those that are spherically symmetric, out of the 0(N 2 ) modes in 
the system are turned on. These can be thought of as the eigenvalues of (jyf 3 . 

Since N is the effective maximal radius for coupling of the two fermions at different 
sites, it makes sense for the wavefunction in the ground state to decrease from a 
maximum value at a = 1 (neighbouring sites) to zero at a = N (sites at distance 
N). For instance, if the wavefunction if) (a) is such that (^(a)) 2 oc N — a, then the 
average distance between sites is 

J \^(a)\ 2 a(2rrada) N 
W f\ij(a)\ 2 (2nada) 2' 1 } 

which is consistent with having a large average distance between the electrons that 
couple. This form of the wavefunction, if) (a) oc y/N — a, here just a consistent choice, 
is exactly what we obtain in the ABJM model. Of course, in principle, if we would 
be able to correctly describe the interactions between various tpf?, as in the ABJM 
model, the dynamics would select the form of ip(a). Finally, the Hubbard model field 
bi must be the linear combination of the spherical modes, i.e. b{ ~ J2 a ^(o)0" a . 

We have already seen that to obtain the Landau- Ginzburg model from ABJM, we 
have only one field, X2, turned on corresponding to turning on the matrix G 2 = 
\/N — m5 m+ i i?l , with (G 2 G\) mn = (N — m)5 mn . As in the simple model above, there 
are two independent rotations, in this case U(N) x U(N) rotations, acting on the 
indices, so the most general solution for the matrix G 2 is in fact U (V ' N — mb mr ^)V~ x . 
We can use these to diagonalize the matrix, thus reducing the degrees of freedom 
turned on, from 0(N 2 ) to 0(N), as in the above condensed matter model. The 
ABJM field that is turned on is X2(G 2 ) mn , corresponding to bi> ~ ^^a^. 

While this field is the only one turned on in our simple toy condensed matter model, 
there are, in principle, many more fields. We could, for instance, have more free 
electrons at each site, thus having more matrix scalars, transforming in some R- 
symmetry group (in ABJM we have 4 complex scalars, corresponding to 0i, 2 , Xi> X2, 
that transform under the SU(2) x SU(2) of the mass deformed ABJM). Then , we 
could also have matrix fermions, corresponding for instance to two electrons at site i 
coupling with one electron at site j, although such modes are, of course, not turned 
on in the Hubbard model description. To complete the field content of the ABJM 
model we need also the Chern-Simons gauge fields, but since those are topological 
and have no dynamics, we don't need to introduce any new degrees of freedom. 

Chern-Simons gauge fields are, of course, no strangers to condensed matter systems, 
showing up, for instance, in the fractional quantum Hall effect (see for instance the 
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review [22]). An abelian Chern-Simons field can be obtained as follows. Consider a 
multi-electron wavefunction ^ e (fi, ...,fk) with a generic Hamiltonian 



g. = E a ~ i ^ w)|! +E''«-3> (5-5) 

3 i<3 

such that H e ^ e = E^ e . We can redefine the wavefunction through the transforma- 
tion 



n 



e -i(pa(ri-fj) 



i<j 



^ e (fi,...,r fc ) (5.6) 



$(ri,...,r fc ) = t/^ e (fi,...,f fc ) 

where a(r*j — r}) is the angle made by = f*j — fj with a fixed axis. Since 

f/- 1 ® - ei?(fi))Z7 —Pi~ eA(fi) - ea(f), (5.7) 

where 

ea(fi) = 0^Vja(fi -r}), (5.8) 

the new Hamiltonian reads 

H = E - — ^ — — + E - ^ ( 5 - 9 ) 

3 i<3 

so that iJ$ = EQ. Therefore after the transformation, a(r) describes a gauge field 
with no dynamics which, one can show is of Chern-Simons type. Such a Chern- 
Simons gauge field, coupled to the fermions and to the electromagnetic gauge field, 
plays a central role in the fractional quantum Hall effect, see e.g. [23]. 

A generalization of this construction to the nonabelian case is straightforward. If 
two fermions at sites i and i" couple to form a boson 0°; a ', at site i' at the midpoint, 
and two other fermions at sites j and j" couple to form a boson 0^' at site j' at their 
midpoint, we can consider the field 

ea{f v ) = Vi> 22 a _ ^i) ' ( 5 - 10 ) 

where we have not yet specified the nonabelian indices on the gauge field. It is not 
hard to see that the only variable in this object is the vector r^/ — fjf (by changing 
the vector we just produce a harmless global spatial translation in the value of 
the right hand side of (5.10)), as well as the discrete choice of fw to belong to the 
fixed point %' or the summed point j'. Since the two vectors f subtracted correspond 
to matrix indices (aa') and (&&'), we can think of this construction as giving us two 
nonabelian gauge fields a ah and a a ' b \ like the A and A of ABJM. Moreover, the 
scalars 0" a ' are bifundamental with respect to the two resulting gauge fields. There 



- 26 - 



remain many open problems to understand about this model, not the least of which 
is the symmetry group acting on the matrix Chern-Simons fields but we leave these 
to the interested reader. This concludes our description of the field content of ABJM 
and qualitative undestanding of abelianization. Suffice it to say that the ABJM 
abelianization gives a well motivated model of AdS/CMT. 

Finally, a few comments on a four dimensional picture for the Landau-Ginzburg 
model (5.1). The Landau-Ginzburg model makes more sense from a theoretical 
viewpoint as a reduction of the corresponding four dimensional theory. But here as 
well, the abelianization presented has in particular an ansatz with the scalar VEV b 
multiplying the matrix G . If we had the same VEV multiplying both G 1 and G 2 , 
that would lead to a description of the fuzzy two-sphere, a finite N approximation 
of the clasical two-sphere [14, 15]. As it is, we can think of the abelianization as 
generating a single direction, or a "fuzzy circle", therefore the resulting Landau- 
Ginzburg theory must also be thought of as coming from a circle reduction of a 
similar theory in 4 dimensions. The physical radius obtained from a fuzzy space 
construction was argued to be (see for instance [14]) 



2 2 
R P h = ^ Tr 



X T X\ =^Tr C J C\ 4vr 2 /f>, (5.11) 



where Ip = l 2 Rn- Assuming this same formula holds for the less defined "fuzzy 
circle" case, from TrfG 1 ^?}] = N(N - l)/2, we get 

^ = (iV-l)|6| 2 47r 2 ^n. (5.12) 

If, as in the pure fuzzy sphere case, the 11th direction has radius Ru = R p h/k, we 
obtain in the "maximally Higgs" case \b\ 2 = c 2 

R ph = (N-l)fil 2 . (5.13) 



6. Some BPS solutions with spacetime interpretation 

We now return to the more general abelianization ansatz, and consider the system 
with 0i = 02 = and gauge fields put to zero, but Xi 7^ still a field (unlike in the 
abelian Higgs case previously described). This ansatz gives the reduced action 

s = - N{N ~ 1) [d*x[\d,xi\ 2 + \dM 2 



2 

4tt 2 



+ ^(lXi| 2 (lx 2 | 2 -c 2 ) 2 + |x 2 | 2 (|x 1 | 2 -c 2 ) 2 ) 
which we now proceed to study. 



(6.1) 
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6.1 Single Profile Solution 



As a first pass, let's consider solutions with a single profile 

Xi = X2 = fW , (6.2) 

with x\ as one of the spatial coordinates. The equation of motion for the reduced 
action (6.1) for this ansatz become 

from which we distill two cases: 

1. Zero mass: In the massless case, \i — 0, the ground state solution is simply 
f{x\) = with no other constant solutions. This is, however, not the only 
solution and a straightforward integration of the equation of motion yeilds 



This is a fuzzy funnel solution which we can check is, in fact, BPS. Indeed, the 
energy of solutions satisfying the above ansatz is 



H^q = -N(N - 1) J dx x dx 2 



2 



(6.5) 



k 2 

which, by the usual procedure of completion of squares, can be expressed as 



H^ = —N(N - 1) J dx 1 dx 2 Ud x J - y/ 3 ^) 



+ surface term I , (6.6) 



from which we can simply read off the BPS equation 

= y/(*0 8 . (6-7) 

It is clear that this equation is solved by the fuzzy funnel solution above. 

2. Nonzero mass: In this case, the constant solutions to the equations of motion 
are 

f = 0, f=M, f=M. (6.8) 

V Z7T V 67T 

Of these, only the first two are ground states. Indeed, completing squares 
again, we find the BPS equation 

d x J^f(f-f]=^ (6.9) 



k V 2n J 
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from which see that indeed / = is a trivial ground state, while the second so- 
lution (f 2 = fik/2n) is again the fuzzy sphere ground state. The third solution 
of the equations of motion (f 2 = fik/6n) doesn't satisfy the BPS equation, so 
is a non-ground state fuzzy sphere. The BPS equation has nontrivial solutions 



fik/2ll 



(6.10) 



/_(0) = +oo, /_(+oo) 



1 =p e - 2 ^i 

The first solution, /_, describes a fuzzy funnel with x\ 6 (0, +oo), so /_ 
varies between an infinite size at x\ — and the fuzzy sphere ground state 
at x\ — > +oo, 

' fik 
2^ ' 

The second solution, /+, describes a fuzzy funnel with x\ G (— oo, +oo) 
in size between zero at x\ — > — oo and the fuzzy sphere at x\ — > +oo, 

' ilk 

This fuzzy funnel solution will be elaborated on in the next section, where we 
argue that it is a generalization of the Basu-Harvey solution that describes an 
M2 ending on a spherical M5. These solutions are plotted in figure 4. below. 



/+(-oo) = 0, /+(+oo) 



(6.11) 
varying 

(6.12) 



/-(si) 



Figure 4: The normalized single-profile solutions f-{x\) and /+(xi) 



6.2 Two-Profile Solution 



The above single profile solution is also fairly easily generalized to a two-profile one 
with 

Xi = f(xi), X2=g(xi). (6.13) 
With this ansatz, the equations of motion reduce to 

2 ~ / ~ jj>k s 



4vr 2 
d 2 f = —f 
1 k 2 



4vr 2 

d ^9 = 



9 2 ~^y -2,y 3 (/ 2 



f 2 -^) 2 + 2f(g 2 



2tt 

fik 
2^ 



(6.14) 
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Again, we can complete squares in the Hamiltonian and read off the BPS equations 



9xJ + f f ( 9 



2 _ (±h\ 

2ir) 



(6.15) 



As in the single profile case above, there are again two separate cases that need to 
be solved separately. 

1. Massless case: For \i = 0, we have the solutions 



k 




Ce 



2vr y/\ - C 2 e~ 2 ^ 



(6.16) 



k 




1 



2vr Vl - C 2 e~ 



2xi 



that solve both the first order BPS equations of motion as well as the general 
second order equations. This solution blows up at x = logC, and goes to a 
constant in g and zero in /, corresponding to a fuzzy circle. These solutions 
are plotted in figure 5. below. 



ii 

^ 0.6 



Figure 5: The normalized two-profile solutions f{x\) and f(x{) 



2. Nonzero mass: For // ^ 0, we have the solutions (see figure 6) 

^Ltfc Cexp [iixx - \e 2 ^ 1 ] 



g[xi 



2vr ,/l^C 2 



exp(— e 



2/ixi 



(6.17) 




2vr ^1 -C 2 exp(-e 2 ^i) 
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Figure 6: The normalized two-profile solutions f(x\) and g{x\) for fj, ^ 0. 

7. Funnel solutions as M2-M5 brane systems 

In this section we will try to find a spacetime interpretation for the fuzzy funnel 
solutions in eq. (6.12). The fuzzy funnel solution (6.4), interpolating between a 
sphere of infinite size and a sphere of zero size, is known to have the spacetime 
interpretation of a flat M2-brane ending on a fiat M5-brane. From the point of view 
of the M2-brane theory given by the massless ABJM, the M5-brane appears as a 
spherical funnel solution, a M5-brane that grows from zero size at x\ = oo to infinite 
size at x\ = 0. We will review this case, reduced to string theory, i.e. a D2-brane 
ending on a D4-brane, later. Also, from the point of view of the M5-brane theory, 
we can write a Blon type solution, corresponding to an M2-brane growing out of the 
M5-brane (directions and 5 are trivial, and in directions 1-4 the M2-brane appears 
as a Blon). From the point of view of the spacetime theory, we have an M2-M5 
system preserving 1/4 supersymmetry in flat space, and we also have an M5-brane 
solution in the (backreacted) background of M2-branes. This picture matches nicely 
with the two worldvolume descriptions. 

With this in mind, we expect that the fuzzy funnel solution of (6.12) should have a 
similar interpretation. The solution interpolating between zero and a fuzzy sphere 
vacuum was found in [24, 16], and we would guess that it can only match with a 
spacetime solution corresponding to an M2-brane ending on an M5-brane. We will 
see however that there is some ambiguity, related to the existence of two solutions, 
the one from zero to the fuzzy sphere, and one from the fuzzy sphere to infinity. 

7.1 Massless case: A Fuzzy Funnel Review 

The solution (6.4) corresponds in spacetime to a flat M2-brane ending on a flat M5- 
brane, a solution which preserves 1/4 supersymmetry as follows. In a flat background, 
the leleven dimensional gravitino transformation law, 

= D^e + W pa \ ~ S^r^iVA , (7.1) 
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must be set to zero in order to obtain a BPS solution. The M2-brane solution ex- 
tended in the (0,l,2)-directions corresponds to a nonzero 3-form A i 2 , with a nonzero 
field strength component F^r (here r is the radial part of all the coordinates trans- 
verse to the M2). The solution is given by a local supersymmetry parameter e(r) 
which is a scalar function of r times a constant susy parameter e( ) satisfying 

r 012 e (0) = ±e {0 ) . (7.2) 

The M5-brane solution extended in the (0,1,3,4,5,6) directions similarly gives a 
nonzero field strength Fg 1 ,,£ 4 , where #i,...,#4 are the four angles obtained for the 
transverse directions (2, 7, 8, 9, 10). Again, the solution for the local supersymmetry 
parameter e(r) is a function of r times a constant susy parameter e( ) satisfying 

1 e(o) - ±e(o) • {<■<>) 

We can then have a solution for an M2-brane ending on an M5-brane preserving 
1/4 supersymmetry by imposing both conditions (which are now compatible). We 
can then reduce this system to 10-dimensional string theory, thereby considering a 
D2-brane ending on a D4-brane. 

From the point of view of the D4-brane theory in a flat background spacetime, the 
fuzzy funnel solution looks like a Blon-type solution. For a spacetime D2-brane in 
the (0,l,2)-directions, called t,x,z, and a D4-brane in the (0,l,3,4,5)-directions, with 
polar coordinates r,9,(f> for the directions (3,4,5), the worldvolume gauge field flux 
on the D4-brane is 

F = (2-Ka')nsm6d6d(f) . (7.4) 

Because the solution is of the Blon type, with the D2-brane growing out of the D4- 
brane, we consider z = z(r) on the worldvolume, leading to DBI D4-brane Lagrangian 

£ = T 4V / (l + ^(r) 2 )(r 4 + (An) 2 ), (7.5) 

where A = 2na'. Since C is independent of z, it follows that dC/dz' is a constant, 
which we can put equal to An, in which case we obtain 

z = ±^ z = — (7-6) 

This corresponds to a funnel solution for a semi-infinite D2-brane ending on the 
D4-brane. 

A similar story takes place for the case where there is a background created by other 
D2-branes (parallel with the first). It is however easier to describe what happens 
in the case of the type IIB solution for D3-branes ending on D5-branes (instead of 
D2-D4), since in that case the spacetime background is easier (there is no M theory 
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reduction). Consider the background generated by other D3-branes, with harmonic 
function /(r), 

ds 2 = f(r)- 1/2 (-dt 2 + dx 2 + dy 2 + dz 2 ) + f{r) +1/2 {dr 2 + r 2 dVt\ + ds 2 ) 

(7.7) 

C( 4) = (f^ 1 - l)dt A dx A dy A dz 

The DBI Lagrangean for the D5-brane reads 

C = T 5 U{1 + f'{r)z'{r) 2 ){r 4 + f-\r)X 2 n 2 ) - An^r)" 1 - l)z'(r)l , (7.8) 

and the same calculation leads to the same solution z(r) = Xn/r, with the function 
f(r) dropping out completely. One can also take the near- horizon limit and consider 
the usual scaling r = a'U, leading to a finite funnel solution that can be interpreted 
from the point of view of the D3-brane theory as 

!/(,) = £. (7.9) 

Note that if we consider a spherical D5-brane ansatz, oriented in the (t, x, y, z, re- 
directions, we obtain the Lagrangian 

C = T 5y /f-i(r)(r* + g-^Xn) 2 ) - An (f-\r) - l) . (7.10) 

If we take the full harmonic function f(r) = 1 + ^, there is no fixed sphere solution 
with r = R =constant, but if we drop the 1 in /, i.e. at very large r, we obtain an 
identity by varying with respect to r, namely Xn/Q — Xn/Q = 0. Therefore in flat 
space, we have asymptotically a solution for very large radius sphere, but at small 
radius we only have the funnel solution; a fixed sphere is not a solution. 

7.2 Massive case: Supersymmetry and a Fluctuation Solution on the M5- 
brane 

The mass deformation changes the 11 dimensional background spacetime from flat 
to [17, 19] 

ds 2 = H~ 2/3 (-dt 2 + dx\ + dx\) + H 1/3 (dx 2 3 + ...dxl ) 

(7.11) 

F A = 2fi(dx 3 A dx 4 A dx 5 A dx 6 + dx 7 A dx 8 A dx 9 A dx 10 ) + dx° A dx 1 A dx 2 A dH' 1 

where H{r) = 1 — |/i 2 r 2 . A naive guess is that the M5-brane has to live in the 
(0, 1, 2, 9, 0, £)- directions, where 9,<f) and£ are the angular directions of (3,4,5,6), 
with r their radial direction, giving 

r 012e< %) = ±6 (0) , (7.12) 
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so that the transverse M2-brane would have to be in the (0, 1, redirections, giving 



r r e (0) = ±e (0) . 



(7.13) 



However we observe that the 4-form in (7.12) has nonzero Fqi2 T as wanted, but since 
Four = 9 r (Aoi2) and not ^(Aoir), there must be a nontrivial Maxwell transformation 
that brings the gauge field A into the desired form. 

How would this M2-M5-brane solution look from the point of view of the D4-brane 
theory (i.e., reducing to lOd string theory and focusing on the worlvolume theory)? 

From the fuzzy S* 2 picture in the ABJM theory an action was found for the fluctu- 
ation modes around the ground state [14]. For the scalar $ corresponding to the 
fluctuations of the radius of the S 2 (transverse direction) the action reduces to just 
a massive mode, i.e. with potential 



Later, in [19] the same fluctuation action was found from the DBI action of a D4- 
brane in the background (7.12). 

For such a potential, a solution was found in [25], and was called the BIGGon (in 
analogy with the Blon), representing, in spacetime, an S 3 giant graviton in type IIB 
on the maximally supersymmetric pp wave background with F-string spikes attached 
at the poles. Since the pp-wave background is T-dual to (7.12) (see, for example, 
[17] for the explicit construction), the same solution should apply in our case. The 
BIGGon was found by similarly taking a single scalar field $ on the 3-sphere (a 
fluctuation of the radial coordinate), with the same potential, but on S 3 instead of 



14 



-[$ 2 + (V 52 $) 2 ] , 



(7.14) 



S 2 , i.e. 



i[$ 2 + (V 5 3$) 2 ] . 



(7.15) 



The solution is 



Q 

sin 




(7.16) 



giving the full radial coordinate (background plus fluctuation) 




(7.17) 



where we have taken the S 3 parametrization to be 



X 4 = R cos ijj , 
X 3 = R sin if; cos 9 , 
X 2 = R sin ip sin 9 sin (f> , 
X 1 = R sin ip sin 9 cos . 
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Therefore in our case we have the BIGGon solution 

= (7.18) 

where the parametrization of the S* 2 is 

X 3 = Rsin9, 

X 4 = R cos 9 sin (ft , 

X 5 = R cos 9cos(f). 

This solution indeed corresponds with our naive expectation of a D2-brane extending 
out perpendicularly from the spherical D4-brane. But the action that it extremizes 
corresponds to small fluctuations of the field $. However, in [19] it was shown how 
to write the full DBI action for D4-branes in the mass-deformed spacetime. 

Massive case: full funnel solution 

In the background (7.12) it was shown that the DBI action for the D4-brane has 
a fixed sphere solution, corresponding to the fuzzy sphere solution of the massive 
AB JM. Now we want to see if we can also find funnel solutions corresponding to the 
ABJM solutions (6.12) and extending the perturbative BIGGon solution above. 

To this end, we consider again an M2-brane in the (0,l,2)-directions and an M5- 
brane in the (0,1, 3,4,5, 6)-directions, with (3,4,5,6) in polar coordinates r,9,<p, and 
£. We reduce M-theory to type IIA on £ and look for a D4-brane extending along 
0, 1, r, 9, <p, with z = z(r, 9, 0) in order to have a Blon-type solution corresponding to 
a perpendicular D2-brane as above. 

Dimensionally reducing the background (7.12) to type IIA string theory, and writing 
only the terms in directions parallel to the D4-brane, we find [19] 

U f H ^ 4- 1 \ 

C (5) = -T^T ( ^ ) dt A Z ' dr Ar4rffi 2 + ••• 

C (3) = (H~ l - l)dt A dx A z'dr + (7.19) 



4 3/2 



#1/4 



kR* 

and a worldvolume flux F = 2XNd^2, giving 



7 = \F - B = 2XN ^r 4 . (7.20) 

2kR* V ; 



Substituting this ansatz into the action 



S = T 4 



d 5 xe-*^-det(g + T) + / (C (5) + C (3) A F) 



(7-21) 
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gives 



$ ■ I ! / df A r/.r A <lr A r/0-> . ( I :'-7/"-') I l + tf- 1 ^-^ ( 2XN '"" 




r 6 V 2fciL 



+ J dtAdxAdrA dn 2 ^XNz'iH- 1 - 1) - ^^z' H ^ + V j j . (7.22) 

Evidently, the Lagrangian £ is independent of which means that dL/dz' is a 
constant, which we can put equal to — T 4 2AiV, in which case 



2AiV _ I M 3 r 3 
r 3 2 ^ 16 



1 4 44 + 42 8 



j /i 2 r 2 

(7.23) 



Here we need to have z' < 0, since the equation is obtained by squaring an equation 
whose left hand side is linear in z' and has a positive coefficient, and whose right hand 
side is negative, after which we take the square root of z' 2 . Since Rq = 2XN[j,kR*, 
after defining 

x = fi 2 r 2 ; y = ftz; a = fiRo, (7.24) 

we obtain the equation 

/ a 2 1 1 x \ f\ x 

(7.25) 



dy \x 2 2 + 4 2 j V 1 4 

dx I -1 x xfl_ 1 x^_ a?_ 

V 1 4 4 4 + 4 2 42 

However, as explained in [19], we are in the approximation a = [iRq <C 1, and the 
fixed sphere ground state solution is r = Rq, or 1 = a 2 < 1. That means that we 
can assume x small in the above equation, thus 



dy a 1 

— ~ h - 

dx 2x 2 4 

which gives finally 

R 2 ^r 2 



2/ir 2 4 

This z(r) has a minimum at 

f4XNkRA 1/4 f y/2R 



UXNkRA 



\ I 1 



(7.26) 



(7.27) 



1/2 

> Rq, (7.28) 



and at r*, the minimum value of z is 



y/fiXNkR* = ?j= . (7.29) 
v2 
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Note that here R is written in physical spacetime variables, R 2 h = 8ir 2 Nlpf 2 , where 
/ = a/ /i/c/(27r) is the radius in the M2-brane worldvolume theory. 

To compare with the fuzzy funnel solutions (6.12), we note that z is the equivalent 
of the M2-brane worldvolume direction x\ and r is the equivalent of the transverse 
direction f(xi). Therefore we must consider r(z) but, since it has two branches, we 
must choose only one. The two branches are: r(z) going from to r* (for z going 
from oo to z m i n ), and r(z) going from r* to infinity. That would naively match the 
two solutions in (6.12), except for the fact that r* ^> R Q , and r* — >■ oo at fi — > 
(with N very large), whereas i? — 2/iAiV — > as \i — > 0. In fact, for // — > 0, we 
obtain from (7.27) 

r « = ^- (7 - 30 » 

This can be compared with the fuzzy funnel solution (6.4), written as 

so the spacetime solution is a deformation of the fi = case. It also matches with the 
first solution in (6.12) in the \i — > limit. However if \x is fixed, the solution becomes 
(7.31) in the X\ — > limit, corresponding to z — > 0, which is not even reachable by 
(7.27). 

It is therefore unclear to us how to relate the two branches of (7.27) to the two 
solutions of (6.12) precisely, other than through the general qualitative behaviour. 
We will leave a precise understanding of the matching to future work. 



8. Conclusions 



In this paper, we have studied various ansatze for abelian reductions of the ABJM 
model, in the general case of nonzero mass, and used them to build a better defined 
AdS/CMT model. We have found a general abelianization ansatz (2.22, 2.21), using 
the matrices G a that describe the fuzzy funnel BPS state and fuzzy sphere ground 
state, and that represent a consistent truncation. A further consistent truncation 
led to a model with topological vortex BPS solutions, but with \<f>\ — > at both 
r = and r = oo while yet another further consistent truncation led to a relativistic 
Landau- Ginzburg model which, depending on the parameter c 2 = fik/(2ir) and on 
the scalar vev b, extrapolates between between the abelian-Higgs model, and a scalar 
4 theory. 

The second abelianization was used to take steps towards a better defined AdS / CMT 
model, since the ABJM model has a gravity dual, and the abelianization corresponds 
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to the collective dynamics of O(N) out of the 0(N 2 ) fields. We also sketched a simple 
condensed matter model for a solid with free electrons that exhibits the same general 
features as the abelianization and leads to a bosonic Hubbard model, which in the 
continuum limit gives the relativistic Landau- Ginzburg system. It will be interesting 
to see if we can make more the model more concrete and elaborate further on its 
relation to ABJM. If successful, our construction provides, in our opinion, a concrete 
embedding of the AdS / CMT correspondence in string theory. 

In the last two sections, we studied various BPS solutions suggested by the abelian- 
ization, finding some generalizations of known solutions. We tried to find a spacetime 
interpretation for the BPS solutions in (6.12) as M2-M5 systems, with partial suc- 
cess. For small fluctuations we succeeded in matching this with the BIGGon solution 
for an M2 ending on a spherical M5, but for the full system we could only match 
only general qualitative behaviour and not the particular solution. It goes without 
saying that more work is needed to understand these solutions. 
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